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ABSTRACT 

Context. Flow-induced instabilities of magnetic flux tubes are relevant to the storage of magnetic flux in the interiors of stars with 
outer convection zones. The stability of magnetic fields in stellar interiors is of importance to the generation and transport of solar and 
stellar magnetic fields. 

Aims. We consider the effects of material flows on the dynamics of toroidal magnetic flux tubes located close to the base of the solar 
convection zone, initially within the overshoot region. The problem is to find the physical conditions in which magnetic flux can be 
stored for periods comparable to the dynamo amplification time, which is of the order of a few years. 

Methods. We carry out nonlinear numerical simulations to investigate the stability and dynamics of thin flux tubes subject to perpen- 
dicular and longitudinal flows. We compare the simulations with the results of simplified analytical approximations. 
Results. The longitudinal flow instability induced by the aerodynamic drag force is nonlinear in the sense that the growth rate depends 
on the perturbation amplitude. This result is consistent with the predictions of linear theory. Numerical simulations without friction 
show that nonlinear Parker instability can be triggered below the linear threshold of the field strength, when the difference in superadi- 
abaticity along the tube is sufficiently large. A localised downflow acting on a toroidal tube in the overshoot region leads to instability 
depending on the parameters describing the flow, as well as the magnetic field strength. We determined ranges of the flow parameters 
for which a linearly Parker-stable magnetic flux tube is stored in the middle of the overshoot region for a period comparable to the 
dynamo amplification time. 

Conclusions. The longitudinal flow instability driven by frictional interaction of a flux tube with its surroundings is relevant to deter- 
mining the storage time of magnetic flux in the solar overshoot region. The residence time for magnetic flux tubes with 2 x 10 21 Mx in 
the convective overshoot layer is comparable to the dynamo amplification time, provided that the average speed and the duration of 
the downflow do not exceed about 50 m s _1 and 100 days, respectively, and that the lateral extension of the flow is smaller than about 
10°. 
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1. Introduction 

Observations of large solar active regions are indicative of an 
organised subsurface magnetic field along the azimuthal (east- 
west) direction, with opposite polarity orientations (from east to 
west or vice versa) in the northern and southern hemispheres. 
Emerging magnetic structures are in a filamentary state, in the 
form of magnetic flux concentrations (flux tubes) of various 
sizes (e.g., sunspots, pores). Flux tubes rise in the convection 
zone, emerge at the surface, and form bipolar magnetic regions, 
which follow polarity rules (Hale's law) and systematic tilt an- 
gles (Joy's law). 

Theoretical studies indicat e that weak ma gnetic fields are 
transporte d by flux expulsio n dSchiisslerl 1 1 984b and convective 
pumping dTobias et al.ll200 lT) to the lower boundary of the con- 
vection zone, where the toroidal magnetic field is amplified 
by radial and latitudinal velo city shear in the solar tachocline 
(see e.g., ISolanki et al.ll2006b . The stably stratified lower con- 
vection zone, in particular the convective overshoot layer, is 
a likely location for the generation and storage of the large- 
scale azimuthal magnetic flux. Numerical studies simulating a 
layer of horizontal magnetic field in the bottom of the solar 



convection zone indicate that an initially uniform field under- 
lying a field-free layer leads to the formation of m agnetic flu x 
tubes by magnetic Rayleigh-Taylor instability (e.g.. Fan l2001l) . 
Following its formation, a toroidal magnetic flux tube can reach 
a mec hanical equilibrium close to the bottom of the convection 
zone dMoreno-Insertis et alj|1992"|) . The emergence of magnetic 
flux tubes driven by magnetic buoyancy and the properties of 
active regions (low-latitude emergence, tilt angles, proper mo- 
tions of sunspots) require azimuth al flux densities of the order 
of 10 s G in the overshoot region dD'Silva & Choudhuril ll993l: 
iFan et alj[l994t ICaligari et al.lll995HSchussledll996l) . The pos- 
sibilities f or the generation of these flux densities have been re- 
viewed bv lSchtissler & Rempell d2002l).ISchussler & Ferriz-Masl 
d2003l) . and lFerriz-Mas & Steinerl d2007l) . 

The average duration of the solar activity cycle is about 
11 years and the total magnetic flux emerging within one cy- 
cle ranges between orders of 10 24 - 10 25 Mx. The amplifica- 
tion of the toroidal magnetic field to its maximum strength in 
about half an activity cycle raises the following question: how 
can the toroidal flux be stored stably for at least a few years in 
the dynamo amplification region? A related problem concerns 
the feedback of magnetic flux loss on the rate of toroidal flux 
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generation. To more clearly understand magnetic flux generation 
and storage in the convective overshoot region, it is important to 
determine and constrain the effects of flows on the stability and 
dynamics of magnetic flux tubes. Here, we focus on: (a) the non- 
linear development of flow-induced flux tube instabilities, which 
can be dynamically significant in the course of toroidal field am- 
plification; and (b) effects of perpendicular flows on the storage 
of toroidal flux tubes. Consequently, some of the flow properties 
prevailing in the overshoot region can also be constrained, by re- 
quiring that flux tubes with field strengths of up to a few times 
10 4 G are stored in the overshoot region for about a few years. 

W e consider the possibility that the fricti on-induced insta- 
bility dHolzwarth et al.ll2007t rHolzwarthll2008l) leads to flux loss 
from the overshoot layer for field strengths below the Parker in- 
stability limit. We refer to such Parker-stable flux tubes as "sub- 
critical" throughout the paper. We investigate effects of flows on 
magnetic flux tubes in mechanical equilibrium to obtain quan- 
titative estimates, and to answer the following question: how 
strongly do (z) finite-amplitu de perturbations of flux tube s by 
perpendicular flows (see also ISchiissler & Ferriz Masll2007l pa- 
per I) and (if) the frictional interaction of the tube with its sur- 
roundings limit the residence time of flux tubes in the overshoot 
region? 

We approach the problem by considering the nonlinearity of 
friction-induced instability (Sect. |2j, including a discussion of 
the effects of finite perturbations on the magnetic buoyancy in- 
stability (Sect. [231 ), and determine the ranges of flow parameters 
that allow the storage of sub-critical flux tubes subject to radial 
flows (Sect. |3). 

The analyses and simulations were car ried out in the frame- 
work of the thin flux tube approximation dSpruitll 19811) . which 
at present is the only existing approach that can deal with the 
small magnetically induced variations in density, pressure, and 
temperature corresponding to the high plasma /? (> 10 5 ) of the 
deep solar convection zone. 



2. Friction-induced instabilities 

The mechanical equilibrium of a toroidal flux tube in the so- 
lar convection zone requires that the plas ma within the tube ro- 
tates f aster than the surrounding medium. M oreno-Insertis et alj 
(11992b demonstrated that a toroidal flux tube rotating initially at 
the same rate as the external medium can reach this equilibrium 
by developing an internal prograde flow. The speed of this "equi- 
librium" flow depends on the magnetic field strength, depth, and 
latitude. A flux tube subject to non-axisymmetric perturbations 
(varying in azimuth) will develop MHD waves that propagate in 
prograde and retrograde directions. If the internal flow speed (in 
the rest-frame of the external medium) is greater than the phase 
speed of the slowest transversal retrograde wave mode, this wave 
mode will be advected in the pr ograde direction, while being 
amplified dHolzwarth et alj|2007l hereafter paper II). The criti- 
cal flux de nsity for the instab ility is lower than that of the Parker 
instability. iHolzwarthl d2008[ hereafter paper III) extended the 
analysis in paper II to toroidal flux tubes, thus including the ef- 
fect of magnetic curvature force and rotation. 

The linear analyses in papers II and III considered Stokes- 
type friction, which is linear in perpendicular velocity. To test 
the results of the linear stability analysis in the nonlinear context, 
we shall consider the drag force, which is quadratic in velocity. 
This test can be implemented by carrying out fully nonlinear 
numerical simulations of flux tubes subject to non-axisymmetric 
perturbations. 
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Fig. 1. Stability diagram for toroidal flux tubes located in the 
middle of the solar overshoot region. The shaded areas indi- 
cate unstable configurations with the grey scale representing the 
growth rate. Dotted lines outline regions of Parker instability. 

2.1. Linear results 

The hydrodynamic drag force exerted by a steady flow with ve- 
locity v ± perpendicular to a straight cylinder is given by 



Fd - —Co 



P eO"±V± 
7lR t ' 



(1) 



where Co is the drag coefficient, p e o is the density of the external 
medium neighbouring the equilibrium flux tube, and R t is the 
cross-sectional radius of the equilibrium tube. 
The Stokes-type friction force per unit volume, 



-ap e0 \ ± , 



(2) 



where the constant a is the friction coefficient, was considered 
in papers II and III to include the friction-type deceleration in 
the linearised equation of motion. In paper II it was shown that 
the critical field strength for friction-induced instability hardly 
depends on the value of a. Comparing Eq. with Eq. ([TJ, one 
finds that the friction coefficient is of the order 



nRt 



(3) 



where the drag coefficient, Co, is set to be unity, bec ause of the 
cylind rical cross-section and high Reynolds number (Ba tchelorl 
fl967h . 

We consider a toroidal flux tube located in the convective 
overshoot region, parallel to the equatorial plane. For the strat- 
ification of the am bient medium, we use a model convection 
zone developed by ISkalev & Stixl d 1991b . which uses a non- 
local t reatment of convection as described bv lShaviv & Salpeterl 
d!973l) . In the model, the overshoot region extends about 10 4 km 
below the base of the convection zone, which is defined as 
the depth at which the convective energy flux changes its sign, 
at about r = 512 Mm. The thickness of the overshoot layer 
corresponds to about 20% of the local pressure scale height. 
Throughout the paper, the terms bottom, middle, and top (levels) 
of the overshoot region are used to refer to the radial positions at, 
respectively, 2000, 5000, and 8000 km above the lower bound- 
ary of the overshoot region, which is at a radius of r — 502 Mm 
(~ 0.72Ro). 

The linear stability analysis of paper III has provided growth 
rates of the friction-induced instability for toroidal flux tubes. In 
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Fig HI growth rates are shown as a function of the field strength 
(Bo) and latitude (Aq) located in the middle of the overshoot re- 
gion, for a = 1CT 7 s _1 and rigid solar rotatiorfj In paper III, the 
value of a was chosen by determining the average perpendicular 
velocity of the mass elements of the flux tube during the initial 
stages of nonlinear numerical simulations for a flux tube with 
cD 3 x 10 21 Mx. The friction-induced instability sets in for 
Bq > 4 X 10 4 G. The onset of Parker instability is in the interval 
6 x 10 4 < B ^ 10 5 G at low latitudes. 

2.2. Nonlinear simulations 

We carried out numerical simulations of a toroidal flux tube 
using a semi-implicit finite-difference scheme developed by 
iMoreno-Insertisi (Il986l) . tha t was extended to thr ee dimensions 
and spherical geometry by ICaligari et alj (Il995l) . The numeri- 
cal procedure is based on the equations of ideal magnetohy- 
drod ynamics in th e framework of the t hin flux tube approxima- 
tion ( Spruit! 19811) . in the form given bv lFerriz-Mas & Schiisslerl 
d!993lll995 ). In the numerical scheme, the flux tube is described 
by a string of Lagrangian mass elements, which move in three 
dimensions under the effects of various body forces. The equa- 
tion of motion for the material inside the toroidal flux tube, in 
a reference frame rotating with the angular velocity of the tube, 
SI, is written as 



latitude cut at X=30° - mid-overshoot 



Dvj _/ B 2 \ 
Pi ~Dt = \ Pi + to) 



(B V)B 

47T 



+p,[g -flx(flxr)] + 2p ; v ; xfi + F D , 



(4) 



where D/Dt = d/dt + v ■ V is the Lagrangian derivative, the 
subscript ; denotes quantities inside the flux tube (the subscript 
e denotes external quantities in the following). The terms on 
the right hand side of the equation are, respectively, the total 
pressure force (gas and magnetic), magnetic tension force, ef- 
fective gravity (including the centrifugal force), Coriolis force, 
and the hydrodynamic drag force, which is given by Eq. ((T). We 
assume rigid rotation in the solar interior. Although rotational 
velocit y shear can be incorporated into the framework o f our 
model (iFerriz-Mas & Schusslerj[l995t ICaligari et all 1 19951) . we 
disregard it for the sake of clarity, to focus on t he fundamen- 
tal effe cts of radial flows on flux tube dynamics. ICaligari et alj 
(119951) found that external shear does not have a significant ef- 
fect on the stability properties and dynamics of flux tubes in the 
Sun. 

The initial value for the cross-sectional radius of each flux 
tube is taken to be R t = 1000 km for all simulations. For Bq = 
10 5 G, this corresponds to a magnetic flux of about 3 x 10" 1 Mx, 
which is typical of a bipolar magnetic region of moderate size 
on the solar surface. 

The flux tube in mechanical equilibrium is initially perturbed 
by small displacements in three dimensions. The azimuthal de- 
pendencies of the initial perturbations in each coordinate k were 
taken to be in the form 



m=l 



£0 sin 



(5) 



where the amplitudes of all modes are equal. In cases of insta- 
bility, the growth rate of the perturbation was determined by an 
exponential fit to the radial location of the top of the growing 
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Fig. 2. Comparison of the analytical solutions for the growth rate 
of the friction-induced instability (lines) with the numerical sim- 
ulations (symbols) for flux tubes at the middle of the overshoot 
region and 30° latitude. The initial perturbation amplitude (£), 
and the Stokes drag coefficient (a) are indicated in the legend. 



loop. Each fit was limited to the time interval corresponding to 
the initial exponential growth phase of the instability. 

We compare the linear results for a range of friction parame- 
ters a, with the results of nonlinear simulations in Fig.|2j which 
shows the growth rates of the friction-induced instability as a 
function of the magnetic field strength, for 30° latitude in the 
middle of the overshoot region, for chosen values of a. The two 
plateaus where the curves for the linear solutions converge cor- 
respond to Parker-unstable regions according to the linear anal- 
ysis (regions enclosed by dots in Fig. []]). Thus the curve for 
a = 10~ 7 s in Fig. [2] can be seen as a horizontal cut through 
Fig. Q] at Aq = 30°. Each set of numerical simulations were per- 
formed with a fixed value of the total initial perturbation ampli- 
tude in the radial direction, Ar = 2\%\, ranging between 51 km 
and 2539 km. The values of a were chosen such that the curves 
correspond to the range of growth rates found in the simulations. 
The simulations exhibit an overall similarity with the linear re- 
sults in terms of the field strength dependence of the growth rate. 
The linear growth rate of the instability based on a fixed value 
of a corresponds approximately to a certain perturbation ampli- 
tude, Ar, as a function of Bo in the linearly Parker-stable and 
frictionally unstable intermediate regime. 

The coefficient or is a measure of the strength of fric- 
tional coupling of the oscillating flux tube with the surrounding 
medium. For a < 10~ 5 s _1 , the growth rate is proportional to 
a, because the frictional coupling facilitates the amplification of 
perturbations. For a > 10~ 5 s , on the other hand, the growth 
rates begin to decrease with increasing a, because too large fric- 
tion impedes perpendicular movements of the flux tube and thus 
the development of overstability. The a values chosen in Fig. [2] 
correspond to the regime in which friction has a destabilising 
effect. The reason for the correspondance between a and ||| is 
the quadratic dependence of the drag force on the perpendicu- 
lar velocity: a larger initial perturbation leads to a larger per- 
pendicular velocity and thus to a larger a, meaning a stronger 
Stokes-frictional coupling. The nonlinear growth of the instabil- 
ity is shown in a supplemented animation (see the online ap- 
pendix, Fig. ED, for a flux tube with B = 7 x 10 4 G, A = 30°, 
and Ar = 5508 km. 
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Table 1. Growth times for various finite perturbations. 
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Fig. 3. Growth time as a function of the displacement amplitude, 
Z?o = 7 x 10 4 G, Aq = 30°, at the middle of the overshoot region. 
The upper horizontal axis shows the flow velocity amplitude that 
leads to the displacement Ar according to the analytical approx- 
imation discussed in Sect. 13. II for m = 5 and |i>J max -26 m s _1 . 



The variation in growth time as a function of perturbation 
amplitude for Bo = 7 x 10 4 G (see Fig.[2]i is shown in Fig. [3] It 
shows that the friction-induced instability grows faster for larger 
perturbations, following a power law. 

There are two reasons for the nonlinear behaviour of the in- 
stability, i.e., that the growth rate is proportional to the initial 
perturbation amplitude: (;) a is proportional to v±, i.e, the drag 
force is nonlinear, (ii) the spatial variation in the superadiabatic- 
ity, 8 = V - V a d, of the external medium along the flux tube 
becomes increasingly effective with increasing |||, such that the 
excessive magnetic buoyancy at the highest location of the tube 
provides an additional upward acceleration. However, the latter 
"^-effect" is not likely to be the dominant source of nonlinearity, 
because otherwise there would be far poorer agreement between 
the linear and nonlinear results in Fig. [2] This point will be in- 
vestigated further in Sect. |2.3l 

We carried out additional simulations for linearly Parker- 
stable configurations with a perturbation of m = 1, by varying 
the amplitude. The results are given in Table [TJ In all cases but 
one for B = 8 x 10 4 G and A = 60°, friction was taken into ac- 
count. For Bo — 3.5 x 10 4 G, the flux tube is linearly stable (see 
Fig.[TJ, and is also stable to the finite perturbations applied in the 
simulations. For stronger fields the friction-induced instability 
sets in, the growth rate then being proportional to the perturba- 
tion amplitude. For a given latitude and perturbation amplitude, a 
higher Bo leads to a more rapid growth, as predicted by the linear 
approach. Owing to nonlinear effects, the initial m = 1 perturba- 
tion initiates perturbations of higher-order modes, in particular 
if the initial amplitude is large. In the case of unstable flux tubes, 
the eigenmode with the shortest growth time dominates the evo- 
lution and causes an increase in the top position of the flux tube. 
During its early evolution, this rise is approximately exponential 
and can be characterised by the e-folding times given in TableQ] 
In those cases marked with asterisks, the top position of the per- 
turbed flux tube decreases with time. Since each eigenmode de- 
cays on its own individual damping time and modes may interact 
nonlinearly, a unique e-folding timescale of the overall decay of 
the top position is not possible. After all, in these cases we find 
that the top position decreases to about one tenth of its original 
value within 10 3 - 10 4 days. The configuration Bo = 8 x 10 4 G 
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Fig. 4. Radial profile of the superadiabaticity in the outer 
40 per cent of the Sun, according to the solar model used 
(ISkalev & Stixlll991h . The arrows and the corresponding circles 
mark the boundaries of the overshoot region. 



and Ao — 60° is very close to the Parker instability boundary 
(see Fig. [TJ, and the result for a sufficiently large perturbation is 
a rapid growth of an instability induced by the steep (5-gradient 
of the external medium along the tube, which we consider in the 
next section. 



2.3. Effect of stratification 

The most important external quantity in determining the stability 
of a toroidal flux tube in the overshoot region is the superadia- 
baticity, 5. Figure|4]shows the radial profile of log \8\ in the outer 
40% of the solar interior, acc ording to the model under consid- 
eration ( ISkalev & Stix|[l99ll) . From the upper radiative zone to 
the lower convection zone, the superadiabaticity increases from 
negative values and changes its sign at about r = Q.HRq. It in- 
creases by a factor of about 2 x 10 5 through the overshoot layer. 
We assume that a toroidal flux tube is deformed in such a way 
that it extends between the boundaries of the overshoot region. 
The parts of the tube that extend to higher layers (smaller \8\) 
will experience larger buoyancy force and thus be destabilised, 
whereas the parts that extend to deeper layers will be stabilised 
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Fig. 5. The difference in superadiabaticity between the surround- 
ings of the top and bottom parts of perturbed flux tubes, which 
are initially located at the middle overshoot zone with Aq = 30° 
and 7 x 10 4 G, the same initial conditions as in Fig. [3] The drag 
force has been taken into account. The numbers indicate growth 
times of the instability for each simulation. 



owing to sufficiently negative superadiabaticity. These nonlinear 
effects may lead to the formation of a rising flux loop. 

Figure [5] shows the initial superadiabaticity difference, A6, 
between the surroundings of the top and bottom parts of the flux 
tube as a function of the initial perturbation amplitude, Ar, for 
the case presented in Fig. [3l(Sect. 12.2b . We note that the relative 
difference A6/6 is lower than unity (log S ~ -6.0). 

To check whether a variation in 6 is responsible for the non- 
linear dependence of the growth rate on the perturbation am- 
plitude, we carried out simulations with the same initial con- 
ditions as given above (Fig. |5), but without the drag force, to 
avoid a mixing with the frictional instability, thus isolating the 
(5-effect. No instability was found within the considered range 
of perturbation amplitudes. This indicates that the <5-effect is not 
responsible for the nonlinearity of the frictional instability for 
log A6 < -6.0. 

In order to constrain the difference of superadiabaticity re- 
quired for the r)-effect to have a significant role in destabilising 
a flux tube, we carried out numerical simulations at the bot- 
tom of the overshoot layer, where the radial gradient of 6 is 
steeper than in the upper overshoot region. Two sets of simu- 
lations were carried out for latitudes Ao = 60° and Aq = 20° 
and Bo - 12.6 x 10 4 G, without the drag force. The two cases 
correspond to linearly Parker-stable configurations, very close 
to the instability boundary, which is about 2 and 4 kG larger for 
the high- and low-latitude cases, respectively. For both cases, 
we find instability for Ar > 1000 km. Figure [6] shows Ad as a 
function of Ar for the high-latitude case. A rising loop is formed 
for Ar > 1000 km, which corresponds to log \A6\ > -5.8. The 
growth rates are about 60 days for Ar 1076 km and 20 days 
for Ar 1537 km, which are comparable with the rise time in 
the convection zone proper. The rising loop is triggered by non- 
linear effects owing to a sufficiently large magnetic buoyancy 
difference along the flux tube, rather than by overstability of in- 
teracting wave modes. 

To summarise the results obtained in this section, we find that 
a superadiabaticity difference of the order of 10 5 is required for 
a subcritical toroidal flux tube to become unstable in the over- 
shoot layer, due to the (5-effect. 



Fig. 6. Same as Fig. El for B = 12.6 x 10 4 G, Ao = 60°, at the 
bottom of the overshoot region, without the drag force. When 
the initial perturbation amplitude exceeds ~ 10 3 km, the flux 
tube becomes unstable owing to buoyancy differences along the 
tube. 

2.4. The slingshot case 

One may expect that rapidly rising flux loops can be formed by 
localised upflows at the top of the overshoot region. An interest- 
ing question is whether such an eruption of a small part of the 
flux tube can originate close to the top of the overshoot region. 
To test this possibility, we carried out a numerical simulation for 
a flux tube at the top of the overshoot region with Bo = 6 x 10 4 G 
and Aq = 40°, which is a linearly Parker-unstable configuration. 
The perturbation was applied with \£ r \ - 1500 km and m = 15, 
which corresponds to an azimuthal extension of 12° for each 
lifted portion. In the course of the simulation, the wave energy 
in high azimuthal modes is gradually transferred into lower-m 
modes. The tube enters into the convection zone with m = 2 
mode, forming two large-scale loops. This is consistent with the 
prediction of the linear stability analysis. 

Small-scale loops can also originate in the bottom of the 
overshoot region, where the 5-effect (Sect. 12.3b is significant. 
A slingshot effect can occur at the radiative zone boundary, if 
a sufficiently strong localised downflow pushes a small part of 
the tube downward. Subsequently, the submerged part would 
be ejected upwards owing to strong buoyancy. We defined ini- 
tial conditions describing localised downward perturbations of 
Gaussian shape with the azimuthal extension ranging from about 
9° to 18° (10 < m < 20) at the bottom of the overshoot region. 
At the beginning of the simulation, the submerged part of the 
tube rises rapidly from the radiative zone boundary. However, 
because of the drag force, which is proportional to the square 
of the perpendicular velocity, it has already been rapidly decel- 
erated within the overshoot region and its translational kinetic 
energy has been partly transferred into MHD waves propagating 
along the tube. 

3. Displacement of a toroidal flux tube by radial 
flows 

We have so far considered the dynamics of toroidal flux tubes 
based on the assumption of spatially perturbed initial configura- 
tions, without an explicit external driving force. In the following 
sections, we investigate the effects of external flows perpendic- 
ular to the tube axis, which displace the tube owing to the drag 



6 



E. Isik & V. Holzwarth: Flow instabilities of magnetic flux tubes IV 



force, given by Eq. ([T). The purposes are (1) to quantify the ef- 
fects of external flows on the subsequent displacement, and (2) 
to test the possibility of storing a toroidal magnetic flux tube with 
fl> ~ 10 21 - 10 22 Mx in the solar overshoot region for a few years, 
which is comparable to the dynamo amplification time. We ex- 
amine the effects of spatially periodic radial flows in Sect. 13.11 
and of localised downflows in Sect. 13.21 



3.1. Azimuthally periodic flow: linear analysis 

We assume that a toroidal flux tube is deformed by azimuthally 
periodic radial flows to such an extent that the drag force is bal- 
anced by buoyancy and magnetic tension. By solving the lin- 
earised equations of motion for the perturbations of a thin flux 
tube, we can derive an analytical expression relating flux tube 
parameters to flow parameters. This relation allows us to esti- 
mate the conditions in which the deformation becomes so large 
that parts of the flux tube enter radiative and/or convection zones. 
These deformations can destabilise a flux tube in the overshoot 
region, e.g., by the (5-effect. 

3.1.1. Solution procedure 

We adopt t he linearised equations of thin magn etic flux tubes 
as given by iFerriz-Mas & Schiisslej (fl993lll995h and apply the 
drag force exerted by an azimuthally periodic flow of the form 



Fd(4>) = — — exp(wn<£). 

7TK t 



(6) 



The azimuthal wavenumber, m, measures the azimuthal exten- 
sion of the perpendicular flow, itjm. We assume that the flux 
tube reaches a stationary equilibrium state, where the drag force 
is balanced by buoyancy and magnetic tension forces, so that the 
time derivative in the momentum equation vanishes. The com- 
ponents of the resulting linearised equations of motion for per- 
turbations in steady state are given in Appendix lAl Eqs. ( lA.ll i- 
( IA.3b . For an azimuthally periodic flow, the components of the 
displacement as a function of azimuth in cylindrical coordinates 
are 



(7) 



£ R = *Xe| R ■ cos m(f>, 

= -3m|,j • sinm^, 

£. = %e£ : ■ cos m<p, 

where the complex amplitudes £r, £ z are given by Eqs. flA.9| >- 
dA.l It . The phase difference between the azimuthal displace- 
ment £0 and the spherical radial displacement £ r is n/2. The az- 
imuthal perturbation leads to diverging flows around tube crests 
and converging flows around the troughs of the tube, to restore 
hydrostatic equilibrium (see Fig. [7] middle panel). Such a flow 
pattern increases the density deficit in the tube crests, thus it has 
a destabilising effect. The azimuthal perturbation does not have 
a significant effect on the radial position of the tube crests, owing 
to (0 the relatively small ratio £<j/£r and (if) the phase relation 
between £ R and 

In the present analysis, we assume that any displacement of 
mass elements along the tube is negligible before the stationary 
equilibrium is reached. This process involves nonlinear varia- 
tions in density, internal flow speed, and magnetic field strength 
as functions of azimuth and time, which are outside the scope of 
our linear approach. We present the justification of this assump- 
tion and discuss its limitations in AppendixlBl 

We begin the analysis by finding the extreme values of the 
velocity amplitude, that displaces parts of the flux tube to the 
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Fig. 7. Geometry of a flux tube deformed by azimuthally peri- 
odic flows, at three initial locations (solid lines): (from top to 
bottom of the figure) the top, the middle, and the bottom of the 
overshoot region. The latitude is 30° and the field strength is 
6 x 10 4 G. The distance from the middle of the overshoot region 
in kilometres is shown on the left axis, and in units of local pres- 
sure scale height (at r = r m \&) on the right axis. The azimuthal 
wavenumber is chosen such that the troughs or crests of the tube 
partially enter the convection or the radiative zones. The radius 
of the tube is 1000 km. Arrows denote the relative strength of 
the resultant perturbation, which is led by the external flow. 



boundaries of the overshoot layer. Figure [7J shows the spherical 
perturbation % r as a function of <p at three depths in the over- 
shoot region for a flux tube with Bq - 6x 10 4 G, Aq = 30°, and 
R, = 1000 km. The velocity amplitudes corresponding to each of 
th ese three depth s were adopted from the convection zone model 
of ISkalev & Stixl d 1 99 It) . The azimuthal wavenumber were cho- 
sen such that in each case the crests or the troughs of the per- 
turbed tube come very close to the boundaries of the overshoot 
layer. We call this azimuthal wavenumber m a i t and the result- 
ing perturbation the critical perturbation. At the bottom of the 
overshoot region, the flow must be very extended (m < 5) to 
destabilise the flux tube. In higher layers of the overshoot re- 
gion, narrower flows (m < 8 and m < 15) can have a desta- 
bilising effect. The critical azimuthal wavenumber required to 
carry the tube to the boundaries of the overshoot region, m cr j t , 
decreases with increasing depth. The relationship between % r , m, 
and the radial location of the tube is determined by 1) magnetic 
curvature force, which increases with m, and 2) the convective 
velocity, which decreases with depth in the overshoot region. 

3.1.2. Parameter study 

The parameters describing the initial condition of the flux tube 
are the field strength, cross-sectional radius, initial radial posi- 
tion, and latitude. The parameters pertaining to the external flow 
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Fig. 8. Variation in the maximal displacement of flux tubes by 
azimuthally periodic perpendicular flows with four different az- 
imuthal wavenumbers, m, as a function of flow speed (Bo = 
6 x 10 4 G, R, = 1000 km). The vertical axes denote the radial 
distance with respect to the middle of the overshoot region (on 
the right axes in units of the local pressure scale height at the 
mid-overshoot region). Solid lines show the maximal downward 
displacement of the troughs of a tube starting in the bottom of the 
overshoot region. Dashed lines show the maximal displacement 
for the troughs and crests of tubes starting in the middle of the 
overshoot region, and the dash-dotted lines show the displace- 
ment of the crests of tubes starting near the top of the overshoot 
region. For clarity, the effect of upward (downward) flows on 
flux tubes near the bottom (top) of the overshoot region is not 
shown. The shaded regions denote the radiative zone below and 
the convection zone above. 



are the azimuthal wavenumber and the velocity amplitude. In the 
following, we shall evaluate Eq. (0 and Eqs. ( |A.91 l- (|A7TTb for 
given sets of parameters, in order to find the perturbation ampli- 
tude as functions of parameters related to the flow and the flux 
tube. 



Dependence on the flow velocity and the azimuthal wavenum- 
ber. Figure [8] shows the variation in r - r m \&, which is the radial 
distance of the crest of the perturbed tube from the middle of 
the overshoot region, as a function of the velocity amplitude of 
the radial flow for various azimuthal wavenumbers. The min- 
imum flow velocity required to advect the flux tube up to the 
convection zone or down to the radiative zone boundaries cor- 
responds to the intersections of the curves with the boundaries 
of the overshoot layer. For a fixed flow speed and increasing m, 
it becomes more difficult to displace the tube crests from the 
equilibrium configuration, because the magnetic curvature force 
increases. For a given m, faster flows lead to larger perturba- 
tions. In the case for which the tube extends between the bound- 
aries of the overshoot layer, the variation in buoyancy along the 
tube can be significant, because the superadiabaticity changes 
by about three orders of magnitude between the upper and lower 
boundaries (cf. Sect. 12.31 Fig. [?}. At the top of the overshoot 
region, a small-scale deformation (a high m) can drive the flux 
tube into the convection zone proper, whereas for the bottom of 
the overshoot region, coherent downflows on a rather large scale 
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Fig. 9. Displacement of flux tubes by radial external flows that 
are periodic in azimuth, as a function of field strength. The left 
panel is for m — 8, and the right panel for m = 20, to show 
the contrast between the effects of flows with large and small 
wavelengths. The radius, the latitude, the initial depth of the flux 
tubes, the vertical axes, and the line styles in the plots are the 
same as in Fig [8] The tube is subjected to an azimuthally peri- 
odic flow with an amplitude of 10, 26, 31 m s _1 for the bottom, 
middle, and top of the overshoot region, respectively. The long 
dashed lines indicate the boundaries of the radiation and convec- 
tion zones. 



(m ~ 5) are required to push the tube down to the radiative zone. 
However, the true depth of penetration close to the radiative zone 
cannot be estimated in this way, because the superadiabaticity 
decreases strongly, and thus the linear approach becomes inap- 
plicable. Indeed, the superadiabaticity becomes strongly nega- 
tive so that the stable stratification largely inhibits any further 
downward penetration (see Sect. I2.41 i. 

Dependence on the magnetic field strength. Figure [9] shows 
the variation in the displacement amplitude as a function of its 
field strength, for R t = 10 3 km, and the amplitude of the az- 
imuthally periodic perpendicular flow is assumed to be 14 m s _1 
at the bottom, 26 m s _1 at the middle, and 31 m s _1 at the top of 
the overshoot region, as in Fig. [7] As the azimuthal wavenumber 
of the flow is increased, the tension force resists the deformation 
of the tube more strongly, so that the perturbation weakens. For 
Bq> 10 5 G and m > 10, the deformation of the tube by the flow 
becomes smaller than 200 km. For m = 20, a perpendicular flow 
with a speed of 26 m s _1 , applied to a tube with Bq = 10 5 G, 
located at the middle of the overshoot region, displaces it by an 
extent of about 5 x 10~ 3 H P (~300 km). The corresponding rela- 
tive change in superadiabaticity between the depths of the crests 
and the troughs is A6 < 10~ 7 (see Fig. [5]), so the ^-effect is inef- 
fective in this case. 



Dependence on the tube radius. We now set the field strength 
to Z?o = 6 x 10 4 G and vary the tube radius, R { . The result- 
ing functional dependence is shown in Fig. [10] for m = 8 and 
m = 20, for tubes starting at the three sets of depths and con- 
vective speeds chosen above for Fig. [9] Because the drag force 
is inversely proportional to R,, thicker tubes are less affected by 
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Fig. 10. Same as Fig. [9] for Bq = 6 x 10 4 G and varying tube 
radius. 



the flow. For m = 20, tubes thicker than about 2000 km in diam- 
eter do not reach the boundaries of the overshoot region. At first 
glance, this indicates that thicker tubes may be stored in the over- 
shoot region for longer times than thinner tubes. However, diam- 
eters larger than about a few thousand kilometres are not relevant 
to the present context, because the thin flux tube approximation 
is not valid if the tube diameter is comparable to the thickness 
of the overshoot layer: the radial variation in superadiabaticity 
across the tube becomes non-negligible, inducing a strong vari- 
ation in the buoyancy over the tube cross-section (see Fig. 0}. 

The upper horizontal axis in Fig. [3]shows the perpendicular 
velocity amplitude of the external flow, leading to a given maxi- 
mum radial displacement of the flux tube, Ar = 2|£, |, for m = 5. 
Using the correspondence between v ± and \£ r \ (Sect. 13.1. IK we 
can make the following estimation, using Fig. [3] A toroidal flux 
tube with A = 30°, B = 7 x 10 4 G, and R t = 10 3 km can be 
stored in the middle of the overshoot region for about 3 years, 
provided that the external flow velocity does not exceed about 
10 m s -1 . This semi-analytical estimate concerning the storage 
of a magnetic flux tube in the overshoot layer is tested using non- 
linear numerical simulations involving a localised downflow, in 
Sect. 13^21 

3.2. Localised radial flow 

We next consider the effects of both a localised radial flow, 
which has a finite extension in azimuth and latitude, and a longi- 
tudinal flow. The former describes, e.g., a convective downdraft 
penetrating into the overshoot layer, and the latter is required to 
define the initial mechanical equilibrium state. 

3.2.1. Linear analysis 

Before treating the nonlinear evolution of flux tubes under the 
combined effects of a localised downflow and the longitudinal 
flow, we develop an analytical approach to help us to understand 
the basic physics. We consider a toroidal flux tube subject to a 
localised downflow in the middle of the overshoot region, and 
describe the downward flow speed by a Gaussian function of 
the azimuth and obtain the stationary solution in a way sim- 
ilar to the one in Sect. I3.ll The description of the flow field 



and the stationary solution for the displacement are given in 
Appendix IA.2I Figure QT] shows the radius-azimuth diagrams in 
which the dashed curves represent the stationary solutions for a 
flow with u max = -10 ms"' (minus sign means that the flow is 
directed downward) and various field strengths. For simplicity, 
the flux tube is located in the equatorial plane. For Z?o = 10 4 G, 
the downflow leads to a valley-shaped deformation in the flux 
tube. In this case, the shape of the flux tube is mostly deter- 
mined by the azimuthal dependence of the perpendicular flow 
speed. For Bo = 3 x 10 4 G and Bo = 4 x 10 4 G, magnetic tension 
increasingly affects the stationary equilibrium shape of the flux 
tube. For a sufficiently strong magnetic field, a localised down- 
flow leads to a deformation with a larger azimuthal extension 
than that of the downflow itself, owing to magnetic tension. 



3.2.2. Nonlinear simulations 

We assume a Gaussian profile for the external flow field, 

~-(r-r m ) 2 (<f>-<p m ) 2 (0-0 m ) 2 



a{t) ■ Umax eX P 



2cr? 



2crl 



2o-2 



(8) 



where a(t) describes the time variation. We set r m = 5.12 x 
10 10 cm and 07 = 5000 km, Q m = 90°, <p m = 180°, so that 
the flow speed has its maximum value, u max , at the upper bound- 
ary of the overshoot region, and falls to about one fifth of its 
maximum in the middle of the overshoot region. The lateral ex- 
tensions and erg are free parameters. For all the simulations, 
the flux tube is located in the middle of the overshoot layer, and 
its cross-sectional radius is 10 3 km. 



Stationary flow (SF) 

To test whether the stationary equilibria predicted in Sect. 13.2.11 
can occur in the nonlinear case, we carried out numerical sim- 
ulations by applying a localised downflow, which becomes sta- 
tionary after a given time, t s = 100 days. We consider a toroidal 
flux tube, which is initially in mechanical equilibrium. We as- 
sume that an external downflow gradually develops around the 
point (r m , 6 m , <p m ), such that the time profile in Eq. ^ is assumed 
to be 



a(f) = 



t/t s 
1 



if t < t s 
if t > t s . 



(9) 



We present two subsets of simulations in this section: for the 
subset SF50, u max = 50 m s _1 , and for SF10 i> max =10 m s _1 . For 
each subset, the free parameters are the latitude, magnetic field 
strength, and the azimuthal and latitudinal extensions (cr^, erg) of 
the downflow, with the condition cr^ = <j g . We define the "rise 
time" of an unstable flux tube as the time between t — and 
the end of the simulation when the tip of the rising loop arrives 
at about r - Q.9%Rq, where the thin flux tube approximation 
breaks down. 

Snapshots from the set SF50 at t — 8 months for various 
initial field strengths are shown with solid curves in Fig. QT| 
in comparison with the linear results for stationary equilibrium 
(Sect. 13.2. It . The direction of the internal flow along the flux 
tube is from left to right. The stationary solution is nearly iden- 
tical to the nonlinear simulation snapshot for Bq = 10 4 G. The 
broadening of the submerged part for higher Bo (see Sect. 13.2. It 
is visible in the numerical solutions. However, the simulations 
deviate from the linear results substantially on the left side of 
the downflow region, to an increasing extent for higher Bq. The 
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Table 2. Rise times for stationary flow simulations. 



Azimuth (degrees) 

Fig. 11. The radial position of an equatorial flux tube as a func- 
tion of the azimuth, as determined by linear analysis (dashed 
curves, u max = -10 ms" 1 ), and by nonlinear simulation of the 
SF50 case, after 8 months (solid curves). The internal flow (in 
the rest frame of the external medium) is from left to right. The 
initial field strength is different in each panel. The initial radial 
position is ro = 5.07 x 10 10 cm (the middle of the overshoot re- 
gion). The dash-dotted horizontal lines indicate the boundaries 
of the overshoot region. 



azimuthally symmetric external downflow deforms the flux tube 
in an asymmetric form, owing to the internal flow. The central 
part of the tube reaches a dynamical equilibrium after about 100 
days, whereas the tube as a whole never comes to equilibrium 
in any case: for low field strengths (1 - 3 x 10 4 G), the up- 
ward portion creates a transverse wave, which propagates left- 
ward (retrograde). In the course of its subsequent evolution, the 
wave energy is transferred to lower azimuthal modes, while the 
submerged portion of the tube persists as long as the downflow 
speed is kept constant. For higher values of Bo, the magnetic 
tension force limits the extent of the downward displacement, 
whereas the excess magnetic buoyancy in the left wing is en- 
hanced. This leads to a Parker-unstable loop for Bq > 5 x 10 4 G, 
because the top of the loop reaches layers of sufficiently high su- 
peradiabaticity, and the density deficit in the loop top increases 
with time. 

To explain the physical reasons for the asymmetric deforma- 
tion, we show in Fig.Q~2]the azimuthal profiles of various physi- 
cal variables at an early stage of the simulation (t - 5 days), for 
the SF50 case, where Bo = 5 x 10 4 G. Once the external down- 
flow begins to advect a part of the tube downward to a smaller 
radius (Fig.[T2h). angular momentum conservation leads to a lo- 
cal increase in the internal flow speed (b). The maximum internal 
flow speed is reached at the left wing of the submerged portion 
of the tube, because of the additional acceleration by the compo- 
nent of gravity along the tube. Through the right side of the sub- 
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merged part, the internal flow decelerates back to its initial equi- 
librium speed. This asymmetric flow leads to a density deficit at 
the left part and a density excess at the right side (c). The re- 
sulting positive buoyancy in the left part increases with time and 
forms an upward moving loop, which is driven further by the 
accelerating internal flow on the right wing of the loop and the 
higher superadiabaticity at the top location (see Fig.[TT] the bot- 
tom panel). Starting from the initial phases, the accelerating flow 
also leads to a slightly lower internal pressure at the left part (e). 
To balance the external pressure, the magnetic field strength in- 
creases on the left part (d), and this contributes to the subsequent 
growth of the magnetic buoyancy instability. For Bq < 5 x 10 5 G, 
the instability does not set in, because the internal flow required 
for the initial mechanical equilibrium is too slow for a sufficient 
density deficit to develop on the left wing of the submerged part. 
For B < 4 x 10 4 G, the transversal wave propagates leftward and 
quits the downflow region, well before a sufficient density deficit 
on the left wing develops. 

Table|2]presents the rise times (in cases of instability) of flux 
tubes at different latitudes, field strengths, and lateral extensions 
of the downflow, for SF50 and SF10 cases. Linearly Parker- 
unstable cases are shown in italics. For the set SF50, the cases 
with field strengths up to 5 x 10 4 G are stable. For the linearly 
Parker-stable cases (B < 9 x 10 4 G), greater lateral extensions of 
the downflow would destabilise the tube, because for larger cr^, 
the magnetic curvature force is smaller, thus the perturbed part of 
the tube can descend deeper, increasing further the internal flow 
speed in the left wing. In this regime, a stronger magnetic field 
also leads to faster growth of the loop, because the internal equi- 
librium flow becomes faster, which develops the density deficit 
on the left wing more rapidly. For the linearly Parker-unstable 
cases (denoted in italics in Table |2), the rise times are gener- 
ally longer, because in most cases the initial central deformation 
of the tube excites standing waves with m — 5. This limits the 
growth of the loop on the left side of the submerged part, and 
the flux tube enters the convection zone only after the growth 
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Fig. 12. Physical quantities illustrating the effect of the external downflow, as a function of azimuth. The snapshots are taken at 
t — 5 days of the SF50 simulation run for fio = 5x 10 4 G. The external downflow is centred around = 180°, which is indicated by 
a vertical line. We plot (a) the radial displacement of the tube from its initial position, (b) the internal flow speed, (c) the difference 
between the internal and external densities, (d) the deviation of field strength from Bo, (e) the internal pressure (dashed line shows 
the external pressure), and (f) the internal density. 



of an oscillatory instability with m = 1, For SF10 simulations, 
the downflow is too weak to trigger nonlinear buoyancy instabil- 
ity, so the linearly Parker-stable cases are stable in the nonlinear 
case, whereas for Bo > 10 5 G Parker instability sets in, in accor- 
dance with the prediction of linear perturbation theory. 

Transient flow (TF) 

Downflows penetrating into the solar overshoot layer are proba- 
bly transient, i.e., they are decelerated by the increasingly stable 
stratification on their way towards the radiative zone. To test the 
consequences of such a flow, we assume a Gaussian profile for 
the time variation in the downflow speed in Eq. ©, such that 



a(t) = exp 



-0 - Uf 



2of 



(10) 



where t m is the time when the maximum flow speed is reached. 
The flow sets in at t — with v = 10~ 3 u max . The initial effect of 
the downflow on the flux tube is similar to that in the SF case. 
However, once the downflow ceases, it leaves behind transversal 
tube waves propagating along the tube, which interfere with each 
other. For unstable cases, the relative longitudinal flow leads to 
instability by frictional coupling, with the growth rate depending 
on the internal flow velocity (thus the field strength) and the per- 
turbation amplitude (cf. Sect. [2J. The difference between SF and 
TF is that in the former case the stationary flow had continued 
for a sufficiently long time to allow the formation of a buoyant 
loop on one side of the downflow region (Fig.fTSI). 

We carried out a systematic survey of numerical simulations 
by varying the flow duration, 2f m , and velocity amplitude, u max 
in the ranges < |u max | < 50 m s _1 and < 2f m < 200 d, for 
Ao = 10°. The constants t m and <r, were chosen such that the 



initial speed is of the order of 10~ 3 u max . For Bo < 5 x 10 4 G, 
the growth time of the friction-induced instability is longer than 
the upper time limit, which we set to be 7 years. Instabilities 
with growth times longer than a few years are not relevant in the 
present context, because within that time the toroidal magnetic 
field must be amplified in the solar tachocline to its full strength. 
For B = 7 x 10 4 G, which is about 20 kG lower than the Parker 
instability threshold for linear perturbations, the flux tube devel- 
ops an unstable loop within about 4 years for some values of 

(|tWl,2f m ). 

The evolution of the tube for \v max \ = 20 m s _1 and downflow 
durations of 60 days (TF60) and 180 days (TF180) is shown in 
supplementary animations (see the online appendix, Figs. IC.2I 
and |C.3l ). The tube rises over about 10 years for TF60, and about 
3.3 years for TF180. In the initial stages of TF180, the evolution 
of the tube resembles the strong-field behaviour for the SF50 
case (cf. Fig.Qj] bottom panel). 

FigureQj] shows the rise times as a function of \v max \ and 2? m . 
Only the final month of the rise time is spent in the convection 
zone proper, the rest being spent in the overshoot region. A rel- 
atively fast downflow (|i> max | < 40 m s _1 ) lasting less than about 
2 months, or a slow downflow (|u max | ^ 20 m s -1 ) lasting less 
than about 6 months allow the storage of the flux tube in the 
overshoot region for more than about 3 years. As a rough esti- 
mate of the flow duration in the overshoot region, we also cal- 
culated the convective turnover time, H p /v m . dx , in the middle of 
the overshoot layer. To lead to the formation of an unstable flux 
loop out of a tube with Bo — 7 x 10 4 G, either the flow should 
last much longer than the turnover time or the flow should reach 
much higher speeds than the range assumed here. 
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Fig. 13. The variation in the rise time for a linearly Parker-stable 
flux tube in the mid-overshoot region, as a function of the du- 
ration and maximum speed of the localised external downflow 
with cr^fi = 5°, 9 m = 180°, r m = 512 Mm (convection zone 
- overshoot layer boundary) for the TF simulations. The field 
strength and latitude of the initial tube are B = 7 x 10 4 G and 
Ao = 10° (cf. Fig.[T}. The diamonds represent the "stable" cases, 
for which no emergence takes place before 7 years, and the size 
of the filled circles is inversely scaled with the rise time. The 
solid curve shows the convective turnover time as an indicator 
for the flow duration, estimated as H p /v m:ix . 

4. Discussion 

We have investigated flow-induced instabilities of toroidal mag- 
netic flux tubes in the solar overshoot region, to extend the re- 
sults of paper III to the nonlinear regime and to quantitatively 
estimate the effects of convective flows on the storage of mag- 
netic flux in the solar overshoot region. Our simulations con- 
firm the main results of the linear stability analysis of paper III 
(Sect. 12. 2\ . The perpendicular velocity component can be ap- 
proximated generally in the form v± ~ 161 ^e(wf), where \£ r \ is 
the amplitude of the radial perturbation and iOf is the eigenfre- 
quency of the fastest growing unstable wave mode. Substituting 
this expression into Eq. (01, we obtain 



C D u ± 
7rR t 



cc |6| fte(w f ), 



(11) 



For a given value of !Re(w t ), which is specified for a given set 
of 6 (depth), Ao, and Bo, a increases with increasing perturba- 
tion amplitude, as shown by numerical simulations in Sect. 12.21 
(Fig. [2}. Considering the relations obtained in Sect. 13.11 a per- 
pendicular flow with u ± =: 1ms" 1 and m = 5 in the middle of 
the overshoot region can displace a flux tube with B — 7 x 10 4 G 
to 61 - 2500 km (cross symbols in Fig.|2|. Substituting u ± into 
Eq. ( fTTT ). this amounts to a ^ 3 x 10~ 7 , which is of the same or- 
der as the value chosen for the full curve in Fig. [2] Considering 
the numerical results in Sect. 12.21 (see Fig. [3}, we conclude that 
an external flow with a velocity amplitude of 1-10 m s _1 perpen- 
dicular to the tube axis leads to instability with a growth time 
of the order of 1000 days (~2.74 years) or l onger. This range 
of velo cities is consistent with the estimates of lvan BallegooijerJ 
d 1982b for the convective velocities in the overshoot region, sup- 
porting the possibility of storing magnetic flux tubes that con- 
tain fluxes of the ord er of 10 21 Mx. In a more recent study, 
iBrummell et al.l d2002l) carried out 3D numerical simulations of 



penetrative convection, which infer penetration depths between 
Q.Q2H p and Q.llH p , when they extrapolate the dimensionless 
numbers in their simulations to solar values. These values are 
also comparable to the radial perturbation amplitudes that we 
have assumed here. 

The numerical experiments presented in Sect. 12.31 have 
shown us that the nonlinear instability occurring in the linearly 
Parker-stable regime is induced mainly by the frictional coupling 
of the flux tube with its surroundings. We have also found un- 
stable flux tube configurations, for which the (5-effect plays a 
significant role in the dynamics. These numerical experiments 
were made without considering the drag force and in the bottom 
of the overshoot layer, where the radial gradient of superadia- 
baticity is steeper than in the remainder of the overshoot region. 
This result may be relevant during the final phases of the de- 
cay of the large-scale toroidal field in the Sun: if we assume 
that the toroidal flux at the upper layers of the overshoot region 
has already been removed to a large extent at this phase, a flux 
tube that forms near the bottom of the overshoot region can be 
destabilised rapidly by strong (possibly rare) convective down- 
flows, leading to a few active regions during a solar minimum. 
We conjecture that the 5-effect is not the main source of flux loss 
from the overshoot layer, with a possible exception in the bottom 
of the overshoot region, provided that overshooting convective 
flows are sufficiently strong. 

In test simulations, we have found that for perturbations 
larger than about 10 3 H p (~55 km), the growth rates start 
to deviate from the predictio ns of linear stability analysis 
dFerriz-Mas & Schiisslerl [19951) . because nonlinear effects gov- 
ern the dynamics and determine the growth rate of instability, 
through either the friction-induced instability (Sect. 12.21 ) or the 
5-effect (Sect. |2~3l ), depending on the depth and the longitudinal 
flow speed. 

In seeking possibilities of a slingshot effect that leads to 
small-scale loops in thin flux tubes, we have found that this effect 
is inefficient in removing magnetic flux from the overshoot re- 
gion. The common result of the simulations is that a loop driven 
by a transient downflow hits the radiative zone and bounces back 
rapidly. However, in its way through the overshoot region it is 
strongly decelerated mainly by friction. 

In a parameter study surveying analytical estimates of 
the displacement amplitude as a function of field strength 
(Sect. 13.1.21 ), we have found that a flow pattern with m = 20 
and u ± = 26 m s acting on a tube with Bo - 10 5 G leads to 
a perturbation of about 300 km. The corresponding difference 
in the superadiabaticity of the stratification is A6 - 10~ 7 (see 
Fig. [5]). Therefore, we do not expect azimuthally periodic flows 
with short wavelengths (m > 20) to trigger flux tube instabilities 
for Bo < 10 G, provided that the displacement is less than about 
300 km, in other words, u ± < 26 m s . 

After understanding the fundamental effects of azimuthally 
periodic flows on flux tubes, we have considered the case of a 
localised downflow in Sect. 13.21 For a given flow pattern, we 
have calculated the evolved states of flux tubes with various field 
strengths corresponding to linearly Parker-stable cases, using the 
steady-state approximation in conjunction with numerical simu- 
lations. The experiments presented in Fig. [TT] can also be inter- 
preted in terms of a toroidal field strength increase with time, 
in the rising phase of solar activity. As the toroidal field is am- 
plified, localised downflows will have stronger disruptive effects 
on flux tubes, owing to the increasing internal flow speed, which 
is determined by the mechanical equilibrium condition. For a 
stationary downflow, we thus suggest that the lower limit to the 
field strength for which flux loops start entering the convection 
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zone is of the order of 5 x 10 4 G for the middle of the overshoot 
region. 

Proceeding to non-stationary, transient flows, we have set 
up a survey of simulations for Bo = 7 x 10 4 G in the middle 
of the overshoot region for Gaussian time profiles (TF case). 
Depending on the dynamical properties of the downflow, we 
have calculated the evolution of a flux tube with an upper time 
limit of 7 years. However, in the first 7 years of an activity 
cycle, the amplification of the large-scale toroidal field in the 
tachocline cannot be neglected. If a flux tube with a sub-critical 
field strength, say Bo = 7 x 10 4 G, remains in the overshoot 
region for a few years, it is likely that it will be amplified and 
eventually form a Parker-unstable loop for Bo ~ 10 5 G, which 
can lead to the formation of active regions with the observed 
properties in the photosphere. 

We have found that the storage of a toroidal magnetic flux 
tube with O 2 X 10 21 Mx for times comparable to the dynamo 
amplification time is possible. Within the scope of the thin flux 
tube approximation and non-local mixing length models of the 
solar interior, we have not found any significant (hydrodynam- 
ically or magnetically induced) stability problem that impedes 
the construction of Parker-unstable tubes with fluxes of the or- 
der of 10 22 Mx. 

We have considered the problem of storing a toroidal flux 
tube field rather than the problem of building up the magnetic 
field in the first place. We have assumed that the field is amplified 
in the overshoot region, e.g., by rotational shear, on a timescale 
of a few years. The stationary equilibria of subcritical flux tubes 
(Bo ^ 10 5 G) advected by radial flows indicate deformations 
comparable to the size of the overshoot region (see Fig. [5). On 
the one hand, it may be argued that these deformations can desta- 
bilise the tube, e.g., by (5-effect, at times comparable to the dura- 
tion of overshooting convective flows (see Sect. 13. It . However, 
nonlinear simulations (Fig.fTTI) show deviations from linear esti- 
mates, as the field strength increases towards 5 x 10 4 G: the tube 
is destabilised as the field strength is increased. These simula- 
tions indicate that toroidal magnetic fields can be stored within 
the overshoot region in the course of their amplification up to the 
critical field strength (10 5 G) required to explain general prop- 
erties of active regions. Our results can be cross-checked when 
the spatio-temporal structure of penetrative convection near the 
bottom of the convection zone is observed by helioseismology, 
or when 3D hydrodynamic simulations with realistic Reynolds 
numbers and sufficiently high resolution become available. 



5. Conclusions 

Based on analytical steady-state approximations and nonlinear 
simulations of toroidal flux tubes in the solar convective over- 
shoot region, we reach the following conclusions: 

- The flow instability driven by the frictional coupling of 
transversal MHD waves with external flows is nonlinear in 
the sense that the growth rate is a function of the initial per- 
turbation amplitude. This is consistent with the results of the 
linear stability analysis of paper III, in which the parame- 
ter a is proportional to the speed of the perpendicular flow 
and thus to the amplitude of the subsequent perturbation. 
Therefore, the perturbation amplitude determines the a pa- 
rameter in the linear analysis. 

- Significant buoyancy variations (5-effect) along magnetic 
flux tubes can lead to a nonlinear buoyancy instability. The 
(5-effect most likely occurs in regions of large radial gra- 
dients of superadiabaticity such as the bottom of the over- 



shoot layer, where a sufficient 6 variation along the tube 
(A6 > 1CT 6 ) is easier to attain than in the upper layers. 
For flux tubes in the solar overshoot layer, we have estab- 
lished links between the perpendicular flow velocity ampli- 
tude, its spatial (azimuthal) extension, and the resulting ra- 
dial displacement and instabilities induced by <5-effect or lon- 
gitudinal flows. 

To store magnetic flux with flux densities between 10 4 -10 5 G 
for times comparable to the dynamo amplification time in the 
convective overshoot layer, the average flow speed and the 
flow duration must not exceed about 50 m s ~' and 100 days, 
respectively, and that the azimuthal extension of the flow is 
lower than about 10°. If we assume that these conditions pre- 
vail in the solar overshoot region, then magnetic fluxes of up 
to 10 22 Mx can be stored within thin flux tubes during the 
dynamo amplification phase. 
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Appendix A: Flux tube subject to radial flows 

The equations for a flux ring at an arbitrary latitude Ao, in the 
limit of B » 1, are given by Ferriz-Mas & Schiissler 1995 
(Sect. 4.1). Consider that a flow along the spherical radial direc- 
tion with an amplitude u ± , is applied to a flux tube with radius 
R t located at latitude Ao- This flow will exert a drag force per- 
pendicular to the tube axis, Fo = Peov 2 ± /nR t , where p e o is the 
unperturbed external density. After a certain time, the drag force 
and the restoring forces of magnetic tension, Coriolis, and buoy- 
ancy will balance each other and a stationary equilibrium will be 
reached. Adopting the equation of mot ion for linearised pertur- 
bations dFerriz-Mas & Schiissler! [l995l) for the stationary force 
equilibrium yields the inhomogeneous system of equations 



d(p 2 



^f + C^ + C R { R + C z { z 

2 f 2 w +D ^ +D ^ +D ^ 



where 



— F D ((f>) cos Ao, 

PeO 



PeO 



-F D ((f>) sin Ao, 



(A.l) 
(A.2) 
(A.3) 



cos Aq\ 2f 
f--^\,A t = --±-smAo, 



A* 



4/ , 
-Ar, Cr — — cos Ao + A cos" Ao, 
J 

IV \ 

C z = 1- A cos Ao sin Ao, 

\ y I 

-- -A z , D R = C z , D z = A sin 2 Ao, 



A = B6 



(A.4) 













2 



(A.5) 

Here we have assumed that the magnitude of the drag force is 
comparable to a first-order perturbation. 
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A. 1 . Azimuthally periodic flow 

We define the flow velocity in the form v± = cos(ra</>) for a given 
azimuthal wavenumber m. The resulting drag force is given by 



F D (cf),m) = Pc ° n ±0 exp(im<p), 
nR, 



(A.6) 



where v ± o is me maximum value of the velocity. The resulting 
displacement is of the form 

% = £ m exp(»n0). (A.7) 
Substituting Eqs. flA.6l & |A!7T i into Eqs. (lA.lb - (lA.3b , we obtain 
A 



' Cr-2 f 2 m 2 imC$ 
im.An —2f 2 m 2 Cj/v 2 ^ 



imA 7 



D R 



imDs 



D z - 2f 2 m 2 ) 



tR 

h 

2 \ 



-^cosAo^ 


v -^-sin/io^ J 



(A.8) 



This equation can be solved in a straightforward way using 
Cramer's rule. The solution for the complex amplitudes yields 



T 2 m 2 v\ Q 



detA nR, 



c 2 

2f 2 ^(D z -2f 2 m 2 )-A z D <f 

vi 



+ sin Aq 



c 2 

A z C^-2f 2 ^C z 

V A 



(A.9) 



XO 



detA nR t 



cosA () [A R (D Z - 2f 2 m 2 ) - A Z D R ] 



+ sin A [A Z (C R - 2f 2 m 2 ) - A R C Z ] , (A. 10) 



r 2 m 2 vj 
detA nR, 



COS Aq 



+ sin Aq 



c 2 

A R D ( „-2f 2 ^DR 



2f^r{C R -2f 2 m 2 )-A R C ] 



(A.11) 



Substituting each component (IA.9b - (IA. lib into Eq. JA.7b . we ob- 
tain 

% R = 2ve(f s>m ) cos(m0), 
^ = -3m(| 0;(n ) sin(m0), 
£ = !Re(|. im ) cos(m0), 



(A. 12) 



A.2. Localised flow 



In the case of a localised perpendicular flow, we define the ve- 
locity and the drag force in the following forms: 



v ± (<f>) = exp 



PeOV 



■ sin 2 (f>/2 



j_o 



nR, 



^]f m exp( 



im<p), 



(A.13) 
(A. 14) 



where g m = g m exp(imcf>) is the Fourier transform of the 
function v^iip). For the resulting displacement, we follow the 
method described in Sect. IA.1I and write the components in the 
Fourier series representation 



(A.15) 



£r = ^ tff,mf m expO'm^), 

m 

£<i> = Zj %</'•">§"> ex p( ! ' m< ^)> 

m 

& = 2 %z,mgm expiirrup), 



where the complex constants ^[R^ Z ] y ,„ are given by Eqs. JA.9b - 
dA. lib - The final solution for the components of the perturbation 
resulting from a localised perpendicular flow has the form 

£ R = ^ 2le(|fi >m ) [3le(| m ) cos(m^) - 3m(g m ) sin(ra0)] , 

m 

= - /^j 3m <&>,>«) [ 3m (lm) cos(m0) + 5Re(£ ffl ) sin(m<£)] , 

m 

4 = 2 ^e(| ; , ffl ) [%e(g m ) cos(m0) - 3m(| m ) sin(m^)] . 

(A. 16) 



Appendix B: The stationary equilibrium 
approximation 

B. 1 . Comparison of azimuthal and radial effects 

Suppose that a toroidal flux tube is in mechanical equilibrium. 
Then apply a localised external flow in the spherical radial di- 
rection. The parts of the tube affected by the flow are advected 
in the direction of the flow, until the drag force is balanced by 
the restoring forces of buoyancy and magnetic tension. To jus- 
tify ignoring any displacement of mass elements along the tube, 
a necessary condition is that the time for the radial displacement, 
t, , is short compared to that of the displacement along the tube 
axis (tangential perturbation, ), so that 



2> 



= TR, 



(B.l) 



where ro is the radial coordinate of the equilibrium tube, ^ is the 
propagation velocity of the perturbation along the tube, and ur 
is the maximum flow speed, and can be determined from the 
solution of the eigenvector problem using the linearised equa- 
tion of motion for first-order perturbations (see Sect. IB. 2b . The 
left-hand side of Eq. ( IB. lb represents a timescale for the re- 
establishment of hydrostatic equilibrium, and the right-hand side 
represents the time it takes for a dynamical balance between the 
drag force and the restoring forces is reached. 

As an example, we consider an equatorial flux tube, in the 
middle of the overshoot layer. The external flow is assumed to 
be localised (Eq. I A. 14b and downwards with a maximum speed 
of 10 m s _1 . This is the same flow configuration as that de- 
scribed in Sects. 13.2.11 and IA.2I For simplicity, we disregard 
the effects of longitudinal wave modes and the presence of az- 
imuthal flow, which is required by the mechanical equilibrium 
condition. Table IB. 11 shows the values for both sides of the in- 
equality ( IB. lb . as a function of the azimuthal wavenumber of the 
external flow, for Bq = 6 x 10 4 G. We also give the values of the 
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downward displacement amplitude, produced by the exter- 
nal downflow. For comparison, numerical simulations were per- 
formed by setting up azimuthally periodic flow configurations 
for various azimuthal wavenumbers. In each case, the flow sets 
in at t — and reaches its full strength (|u max | = 10 m s _1 ) in 
30 days. The tabulated values are the maximum displacement at 
t = 240 days, for which the drag force exerted by the external 
flow is already balanced by buoyancy and magnetic tension. The 
variation in displacement amplitude as a function of m for az- 
imuthally periodic flows shows general agreement with that for 
the linearly estimated amplitudes of each azimuthal wave mode 
corresponding to a localised downflow. 

The ratio t^/tr decreases for smaller m. When becomes 
comparable to tr (as the ratio falls below about 10), the az- 
imuthal perturbation might have a significant effect on the den- 
sity distribution along the tube, such that our assumption of a 
stationary equilibrium between the normal forces becomes inac- 
curate. 

Table B.l. Characteristic times for radial and azimuthal dis- 
placements. 



m 


£r (linear) 
(km) 


£ R (num.sim.) 
(km) 


is(d) 


i>(d) 


Vt r (d) 


3 


14077 


12180 


16.3 


85 


5.2 


4 


4347 


4037 


5.0 


202 


40 


5 


2302 


2366 


2.7 


290 


109 


6 


1461 


1531 


1.7 


358 


211 


7 


1021 


1044 


1.2 


408 


345 


8 


757 


731 


0.9 


445 


509 



B.2. Flux tube in the equatorial plane: eigenvector problem 

For simplicity, we consider a flux tube in the equatorial plane, to 
obtain eigenvectors for the wave-like solutions of the perturbed 
flux tube. In this case, the momentum equation for linear pertur- 
bations is simpler than for a non-equatorial flux tube, and there- 
fore it is simpler to make a physical interpretation. 

To obtain relations between the flow velocity and the result- 
ing displacement, we consider the linearised equations of motion 
for perturbations (£«, where R is the distance from the ro- 

tation axis, cp is the longitude, and A is the latitude. Here, only 
radial and azimuthal components (perturbations in the (R, (f>) 
plane) are considered, because the equation for the latitudi- 
nal component is decoupled in the equatorial plane. Following 
Ferriz Mas & Schiissler (1993), the linearised equations of mo- 
tion for the perturbations in the (R, <f>) plane read 



i*(&-20&) = 2/^-4/1/ 



2y) d(po 



(B.2) 



where 

T = 2(cr-l)/ 2 + i/--(--i 

r r \r 2 
+P5 + t 2 (ct - i)(n e 2 () - Q 2 ), 

t a (Elf -_ V2^, 

VA 



(B.4) 
(B.5) 



/ = 



P 



PiO 



(B.6) 
(B.7) 



The meanings of the various symbols in the above equations are 
as follows: H p is the local pressure scale height, 5 is the supera- 
diabaticity, and p,o are the equilibrium values for the pres- 
sure and density inside the tube, respectively, ro is the spherical 
radial position of the equilibrium tube, O and D. M are, respec- 
tively, the angular velocities of the plasma in the tube and the 
external medium, va is the Alfven velocity, Bo is the magnetic 
field strength, and go the gravitational acceleration with a radial 
dependence goer"", where cr is taken to be -2, which means that 
the small contribution of the mass inside the convection zone is 
neglected. The substitution of the ansatz 



%m = £mexp(;wf + invp). 

into Eqs. (!B.2b -( lB~3l leads to the system of equations 

1 



(B.8) 



(2/ 2 m 2 -T- <h l )i R + i 



4/m / 



2y 



2Q.U) 



= (B.9) 



2Qw 



$ R + (2fW - u 2 % = 0. (B.10) 



For the existence of non-trivial solutions, the characteristic de- 
terminant must vanish. This yields the dispersion relation, which 
is a 4th-order polynomial in terms of the eigenfrequency. The 
eigenvector in the azimuthal direction is determined by one of 
the equations in the system (IB.lOl i: 



2£loj-4fm(f-l/2y) 



(B.ll) 



o> 2 - 2pm 1 

from which we obtain, by using the ansatz Eq. ( IB. 81 ), the velocity 

= iu) m if (B.12) 

of the azimuthal perturbation, as an estimate of the rate at which 
the hydrostatic equilibrium along the perturbed flux tube is es- 
tablished. We use this result in evaluating Eq. dB.lt , to check for 
the validity of the assumption of stationary equilibrium of a flux 
tube being displaced by an external flow. 
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Fig.C.l. The shape of a flux tube with B Q = 7 x 10 4 G, A = Fig.C.3. Same as Fig. IC.2I at t = 180 d, when the transient 
30° and Ar = 5508 km at t — 1030 d, in the (r, <p) plane. The downflow ceases, 
horizontal lines mark the boundaries of the overshoot region. 
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Fig. C.2. The shape of a flux tube with B = 7 x 10 4 G, A = 10°, 
at t = 60 d, when the transient downflow (v max — 20 m s _1 ) 
ceases. The horizontal lines mark the boundaries of the over- 
shoot region and the vertical line denote the azimuthal location 
of the centre of the downflow. 



Appendix C: Animations 

We provide three animated GIF files showing the simulated evo- 
lution of perturbed flux tubes, available on-line. 

The animation frict_inst.gif shows the final phases of 
the development of the friction-induced instability in the over- 
shoot region, between t = 700 d and t = 1100 d (Se ct. l2~2l . 
A snapshot from the animation is shown in Fig. IC.ll The as- 
terisk signs represent selected mass elements, which all move 
rightwards in the initial phases, owing to the internal equilib- 
rium flow. We assume that B = 7 x 10 4 G, A = 30°, and the 
initial radial perturbation amplitude is Ar = 5508 km. The initial 
location is the middle of the overshoot region. 

The animation files TF60.gif and TF180.gif show the 
evolution of a flux tube subject to a radial downflow with a 
duration of 60 and 120 days, respectively, until t = 300 days 
(Sect. [3~2~!2l . Two snapshots are shown in Figs. |C.2| and |C.3l cor- 
responding to the time when the downflow ceases in each case, 
i.e., at t — 60 d and t = 180 d. Note that for the longer-duration 
flow, the flux tube is radially more disturbed when the action of 
the downflow is finished. 



